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ABSTRACT 

Gamma-ray spectra from cosmic-ray proton and electron interactions with dense gas clouds 
have been calculated using a Monte Carlo event simulation code, GEANT4. Such clouds are 
postulated as a possible form of baryonic dark matter in the Universe. The simulation fully tracks 
the cascade and transport processes which are important in a dense medium, and the resulting 
gamma-ray spectra are computed as a function of cloud column- density. These calculations 
are used for predicting the Galactic diffuse gamma-ray spectrum which may be contributed by 
baryonic dark matter; the results are compared with data from the EGRET instrument, and used 
to constrain the fraction of Galactic dark matter which may be in the form of dense gas clouds. 
In agreement with previous authors, we find useful constraints on the fraction of Galactic dark 
matter which may be in the form of low column-density clouds (E < lOgcm"^). However, this 
fraction rises steeply in the region E ^ 10^ gcm^^, and for E > 200 gcm^^ we find that baryonic 
dark matter models are virtually unconstrained by the existing gamma-ray data. 

Subject headings: gamma-rays — cosmic rays — ISM: clouds — dark matter 



1. Introduction 

The nature of dark matter remains one of the 
outstanding questions of modern astrophysics. 
The success of the cold dark matter cosmological 
model (albeit with "dark energy" now required: 
ACDM) argues strongly for a major component of 
the dark matter being in the form of an elementary 
particle. However, the inventory of baryons which 
we can observe locally falls far short of the total in- 
ferred from observations of the cosmic microwave 
background fluctuations (Fukugita 2004) , leaving 
open the possibility that there may be a significant 
baryonic component of dark matter. Furthermore, 
although ACDM is very successful in describing 
the growth of structure in the universe on large 



scales, we still lack a direct detection of any of 
the candidate dark matter particles. Lacking this 
decisive piece of observational evidence, some au- 
thors have proposed models which include a large 
component of baryonic dark matter. In particu- 
lar there have been many papers dealing with the 
possibility that cold, self-gravitating molecular 
clouds constitute a major component of the dark 
matter (Pfenniger, Combes and Martinet 1994; 
De Paolis et al. 1995; Henriksen and Widrow 
1995; Gerhard and Silk 1996; Combes and Pfen- 
niger 1997; Walker 1999; Sciama 2000a,b). A 
variety of different forms, including isolated, clus- 
tered, and fractal, have been considered for the 
clouds, but all proposals involve dense gas of high 
column-density, in contrast to the diffuse gas in 
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the interstellar medium which is easily detected in 
emission and/or absorption. 

One of the fundamental predictions of a model 
featuring dense gas clouds is the gamma-ray emis- 
sion resulting from cosmic-ray interactions within 
the clouds (De Paolis et al. 1995; Kalberla, 
Shchekinov and Dettmar 1999; Wardle and 
Walker 1999; Sciama 2000a). Because of the po- 
tentially large total mass of gas involved, this pro- 
cess may yield a diffuse flux in the Galactic plane 
comparable to the flux from known sources for 
photon energies around 1 GeV (Sciama 2000a). 
Considering the high quality data on diffuse emis- 
sion acquired by the EGRET detector aboard the 
Compton Gamma Ray Observatory (Hunter et 
al. 1997), it is worth considering this source of 
gamma-ray emission in detail as it is possible to 
use these data to constrain the dark matter models 
(see Salati et al. (1996); Gilmore (1994)). Most 
previous investigations of this problem have ne- 
glected the self-shielding and cascade phenomena 
which can be important at high column densi- 
ties (Kalberla, Shchekinov and Dettmar 1999; 
Sciama 2000a), and have employed emissivities 
appropriate to the low-density limit. These effects 
alter the emergent gamma-ray spectrum, and we 
note that this could be relevant to the observed 
excess Galactic flux above 1 GeV (Hunter et al. 
1997). We have noted elsewhere (Walker, Mori 
and Ohishi 2003) that massive (M > 10^ Mq) ag- 
gregates of dense gas clouds could potentially ac- 
count for many of the unidentified discrete sources 
detected by EGRET (Hartman et al 1999). 

Here we present detailed calculations of the 
gamma-ray spectra arising from cosmic-ray inter- 
actions with dense gas clouds. We have used a 
Monte Carlo simulation code, GEANT4, devel- 
oped for simulating interaction events in detec- 
tors used in high-energy particle physics. Not sur- 
prisingly, we find that the predicted spectra dif- 
fer substantially between high and low column- 
density clouds, and we discuss the interpretation 
of our results in the context of the observed Galac- 
tic gamma-ray emission. Our calculations are un- 
dertaken for cold, dense molecular gas in clouds 
of radius R ~ 10^'^cm, similar to those proposed 
by Walker and Wardle (1998) to explain the ex- 
treme scattering events (Fiedler et al. 1994) 
during which compact extragalactic radio sources 
are magnified and demagnified as a plasma "lens" 



moves across the line of sight (see McKee (2001) 
for a criticism of this model). However, the re- 
sults of our calculations depend primarily on the 
column-density of the individual clouds, E, under 
consideration, and their fractional contribution to 
the Galaxy's dark matter halo, and our results can 
be taken as representative of other models which 
are characterised by similar values of these quan- 
tities. 

2. Results 

2.1. Gamma- ray production in dense gas 

Previous calculations of gamma-ray spectra 
from cosmic-ray irradiation assumed single inter- 
actions of protons with the interstellar medium 
(Stephens and Badhwar (1981); Dermer (1986); 
Bertsch et al. (1993) and references therein). In 
order to investigate cosmic-ray interactions with 
dense gas, where cascade processes and particle 
transport are important, we have used a Monte 
Carlo code, GEANT4,^ to derive gamma-ray pro- 
duction spectra. This code is a general purpose 
Monte Carlo code of particle interactions and is 
widely used for simulation of high-energy parti- 
cle detectors in accelerator experiments. Cross- 
sections and interactions of various hadronic pro- 
cesses, i.e., fission, capture, and elastic scatter- 
ing, as well as inelastic final state production, are 
parametrized and extrapolated in high and low 
particle energy limits, respectively. The tt^ pro- 
duction in this code, which is important because 
of the TT^ 27 decay that dominates the emis- 
sivity of the gas at high energies, has been tested 
against accelerator data (Dannheim 1999). 

Initially we experienced one slight difficulty in 
applying GEANT4 to our physical circumstance: 
the low-energy hadron interaction code, called 
GHEISHA, did not conserve energy very accu- 
rately (GEANT4 bug reports No. 171 and 389). 
A "patch" was available for GHEISHA (GEANT 
ver. 4.4.1), but this patch appeared to introduce 
further problems of its own in the energy de- 
position distribution (GEANT4 bug report No. 
415). These difficulties have been overcome by 
the GEANT team, and we are not aware of any 
such problems in the latest release (GEANT ver. 
4.5.1). 

^Available at http://wwwinfo.cern.ch/asd/geant4/geant4.html 
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Our calculations assume a spherical cloud of 
molecular hydrogen of uniform density and tem- 
perature (10 K). The radius of the sphere was as- 
sumed to be i? = 1.5 X 10^"^ cm ~ 1 AU. Protons 
and electrons are injected randomly at a surface 
point of the cloud and particles subsequently em- 
anating from this surface arc counted as products. 
The adopted spectra of cosmic-ray protons and 
electrons were taken from Mori (1997) (here we 
use the "median" flux; note that the units on his 
equation (3) should read cm~^s~^sr~^GeV~^), 
and Skibo and Ramaty (1993), respectively. 
In Figs. 1 and 2 we plot these spectra together 
with observational data. (Note that the adopted 
spectra are those for cosmic rays in the Galaxy, 
whereas the measured points are subject to the 
modulating influence of the magnetic field of the 
solar wind). The simulated range of kinetic energy 
of cosmic rays is from 10 MeV to 10 TeV. We di- 
vided this energy range into four and superposed 
the resulting spectra with appropriate weight fac- 
tors in order to increase the simulation statistics 
at higher energies, considering the rapidly falling 
spectrum of cosmic rays. The density of molec- 
ular hydrogen, g, was varied from 5 x 10""^^ to 
5 X 10~^ g cm~'^ in factors of 10. This corre- 
sponds to the column density, S = 2gR{cos6), of 
10~^, 10^^, . . . 10"'' g cm~^, respectively, where 9 
is the incident angle of a cosmic ray into a cloud 
and (cos^) = 2/3 for random injection. 

2.2. Simulations 

Figure 3 shows the resulting gamma-ray spec- 
tra obtained in our GEANT4 simulation for pro- 
ton injection into a cloud of column density E = 
100 gcm~^. (The quantity plotted is E'^S, with £ 
being the emissivity as given in §3.) The dashed, 
dotted and dot-dashed lines show the spectral 
components classified by the parent processes pro- 
ducing gamma-rays, i.e., 7r° decay, bremsstrahlung 
and positron-electron annihilation, respectively. 
The latter two components are necessarily omit- 
ted in calculations which assume single interac- 
tions only (i.e. the thin material limit). The er- 
ror bars are calculated from Monte Carlo statis- 
tics. Although the 7r° decay component shows a 
broad peak at ~ 70 MeV (note that the quan- 
tity plotted is the emissivity multiplied by E^) 
and dominates above about 200 MeV, the electron 
bremsstrahlung component broadens the emissiv- 



ity peak. The bremsstrahlung and annihilation 
contributions are negligible in the limit of small 
column-density; their fractional contribution to 
the emission is greatest for clouds which are just 
thick enough (S ~ 10^ — 10^ gcm~^) to attenu- 
ate the bulk of the incident proton power. This is 
as expected considering that the interaction mean 
free path of GeV protons is ~ 80 g cm~^ and is 
in accordance with the result of Umebayashi and 
Nakano (1981) who treated a similar problem by 
solving one-dimensional transport equations. 

The resulting gamma-ray emissivities for clouds 
of various column densities are shown in figure 
4 (proton injection) and figure 5 (electron injec- 
tion). Here the emissivities are defined for ir- 
radiation by cosmic-rays of all species (see §2.3, 
equation (3)); to take account of the contribution 
of heavier nuclei than helium, the emissivity due 
to proton irradiation (figure 4) has been multi- 
plied by a nuclear enhancement factor (Cavallo 
and Gould 1971; Stephens and Badhwar 1981; 
Dormer 1986; Gaisser and Shaefer 1992) of 1.52 
(Mori 1997). Note that for high densities the 
Monte Carlo statistics are rather poor, since the 
yield itself is low. Figure 4 includes a comparison 
of our calculated gamma-ray production functions 
with that of Mori (1997) (corresponding to the 
"thin material" limit). The results are consistent 
with those of Mori (1997) for column densities 
less than about 10 g cm~^, except in the energy 
range E > 10^ MeV where the effect of the max- 
imum energy assumed in the Monte Carlo sim- 
ulation is evident. We note the very low values 
of the emissivity at energies > 100 MeV, for col- 
umn densities E > lO^gcm"^. A slightly steeper 
spectrum in the 10^-10^ MeV region comes from 
our omission of the contribution of heavy nuclei, 
which were taken into account in Mori (1997). 
A somewhat surprising feature of these curves is 
that the power-law index above 1 GeV is almost 
the same as the input cosmic-ray proton flux for 
column densities less than about 1000 g cm~^ (for 
higher column densities the statistics of the sim- 
ulations are not good enough to decide whether 
this result still holds). This is already indicated 
by Umebayashi and Nakano (1981), but is con- 
trary to the expectation of Sciama (2000a) who 
suggested a spectral change above 1 GeV at the 
point where self-shielding becomes important. 

For some purposes the energy-integrated emis- 



3 



sivitics arc of more interest than their differential 
counterparts, so we present these quantities in fig- 
ures 6 and 7, for cosmic ray protons and electrons, 
respectively. At low column densities, where cas- 
cades and self-shielding are unimportant, there is 
very little variation of the integrated emissivity 
with cloud column density, and the thin mate- 
rial limit can be adopted for columns less than 
10 gcm~2 for protons (1 gcm"^ for electrons). 
Above this point, however, the emissivity drops 
rapidly with increasing cloud column-density. In 
order to gauge the sensitivity of these calculations 
to the assumed cosmic-ray spectra, we have com- 
puted our results for two different incident cosmic- 
ray proton spectra, and three different incident 
cosmic-ray electron spectra, as shown in figures 
6 and 7, respectively. The differences are seen to 
be small in comparison with the variation as a 
function of cloud column density, but at a fixed 
column-density the systematic uncertainties asso- 
ciated with the input cosmic-ray spectra are nev- 
ertheless significant. 

2.3. Computation of gamma-ray intensi- 
ties 

Using the gamma-ray production spectra ob- 
tained in the previous section, we have calculated 
the diffuse gamma-ray emission from the Galaxy 
as follows. The predicted gamma-ray spectrum for 
each case is 

1 f°° AN 
lD = ^j^ dsp{s)Ms)— (1) 

where dN/dE is the spectriim returned by the 
simulation in units of photons/MeV/primary, for 
an individual cloud, appropriate to the incident 
cosmic-ray spectrum. The quantity Jcr(s) is the 
intensity of cosmic rays at a distance s along the 
line of sight, in units of primaries cm"^ sr~^, 
and p{s) is the mean density in gas clouds of col- 
umn density S. 

The Galactic variation of the spectrum Jcr(s) is 
not well constrained by existing data, and conse- 
quently we adopt the simplifying assumption that 
the shape of the cosmic-ray spectra (both elec- 
trons and protons) is the same everywhere in the 
Galaxy, with variations only in the normalisation. 
With this assumption it is convenient to recast the 
calculation as 

In = ^£Q, (2) 



Table 1: Integrated gamma-ray emissivities of dense 
gas clouds, due to cosmic-ray hadrons and electrons; 
these results assume Mori's (1997) "median" pro- 
ton spectrum and Skibo and Ramaty's (1993) elec- 
tron spectrum. A nuclear enhancement factor of 
1.52 has been applied to the emissivities computed 
for cosmic-ray protons in order to account for heav- 
ier nuclei. Column-densities (S) are given in units 
of g cm~2 and integrated emissivities in units of 
photons s^^g^^; photon energies {E) are given here 
in MeV. 
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Table 2: Gamma-ray emissivities of dense gas clouds, 

due to cosmic-ray hadrons and electrons; these re- 
sults assume Mori's (1997) "median" proton spec- 
trum, and Skibo and Ramaty's (1993) electron spec- 
trum. A nuclear enhancement factor of 1.52 has 
been applied to the emissivities computed for cosmic- 
ray protons, in order to account for heavier nuclei. 
Column-densities (S) are given in units of gcm"2^ 
and emissivities in units of MeVs-^g^^; photon en- 
ergies (E) are given here in MeV. The quantity 
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(phcm ^sr ^MeV where the emissivity is 

47r , ,diV 

(phs~"'^MeV~"'^g~'^) with Jcr(0) the cosmic-ray 
mean intensity in the Solar neighbourhood and 

is the weighted column density (gcm~^) of the 
cloud population along the line-of-sight under con- 
sideration. This formulation is convenient because 
the emissivity, £, describes the properties of the 
gas clouds themselves and is independent of the 
Galactic variations in mean dark matter density 
and cosmic-ray density; conversely the quantity 
Q characterises these properties of the Galaxy, 
and is independent of the properties of the gas 
clouds themselves. The emissivity shown in fig- 
ures 4 and 5 is the quantity E'^S, whereas figures 
6 and 7 show J^dE'E. For the inner Galactic 
disk, where we are interested in {Id), we need to 
average over the whole solid angle, fl, under con- 
sideration: (Q) = J dn Q /n. In order to calculate 
Q we need to adopt models for both the Galactic 
cosmic-ray distribution and the Galactic distribu- 
tion of the clouds. 

The quantity p(.s), the density in cold, dense 
gas clouds, is only weakly constrained by direct 
observation, because the hypothetical clouds con- 
stitute a form of dark matter. We therefore pro- 
ceed by adopting a conventional dark matter den- 
sity distribution for the Galaxy, namely a cored 
isothermal sphere, as our model cloud density dis- 
tribution, with a fiducial normalisation which is 
equivalent to the assumption that all of the dark 
matter is in the form of dense gas clouds. This 
corresponds to the model 
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in terms of cylindrical coordinates {R,z), with 
a = 155 kms"'^. We have adopted a core ra- 
dius of Tc = 6.2 kpc based on the preferred model 
of Walker (1999). (This choice corresponds to 
Walker's preferred value of cloud column density 
E = 140gcm~2.) Walker's model exhibits a core 
radius which is a function of cloud column density, 
but we have fixed the core radius at 6.2 kpc for all 



of our computations. This choice permits more 
straightforward consideration of the observational 
constraints because Q is independent of S in this 

case. 

It then remains to specify the cosmic-ray 
energy-density as a function of position in the 
Galaxy. Webber. Lee and Gupta (1992) (hereafter 
WLG92) constructed numerical models of cosmic- 
ray propagation in the Galaxy; they did not give 
any analytic forms for their model cosmic-ray dis- 
tributions, but an appropriate analytic approx- 
imation can be deduced from the results which 
they obtained. They found that the cosmic-ray 
radial distribution reflects, in large part, the radial 
dependence of cosmic-ray sources, with a modest 
smoothing effect introduced by diffusion. We have 
therefore adopted WLG92's preferred model (their 
model 3) for the radial distribution of sources as 
our model for the radial distribution of cosmic 
rays. 

The various spectra of cosmic-ray isotope ratios 
considered by WLG92 favour models in which the 
diffusion boundaries are in the range 2 — 4 kpc 
above and below the plane of the Galaxy. We 
adopt the midpoint of this range. WLG92 do not 
give a simple functional form for the vertical vari- 
ation of cosmic-ray density within this zone, so 
we have simply assumed an exponential model: 
exp{-\z\/h). We know that in WLG92's mod- 
els the cosmic-ray density is fixed at zero at the 
diffusion boundaries, and consequently the scale 
height of the exponential should be approximately 
half the distance to the diffusion boundary, i.e. 
/i ~ 1.5 kpc. 

These considerations lead us to the model 
cosmic-ray mean intensity distribution Jcr{R, z): 

(6) 

in terms of cylindrical coordinates {R,z). Here 
Ro — 8.5 kpc is the radius of the solar circle, while 
L = 7 kpc, h = 1.5 kpc and Jcr(0) = Jcr(-Ro)O). 
This distribution has the character of a disk with 
a central hole. 

These models for p{s) and Jcr allow us to 
compute the quantity Q, as per equation 4, and 
the resulting variation over the sky is plotted in 
figure 8. For reference we give the values of 
Q{1, b) evaluated at the cardinal points, as follows: 
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Q(0,0) = 6.47X 10-2gcm-2, Q(±90°,0) = 1.54 x 
10-2gcm-2, Q(180°,0) = 7.70X 10-3gcin-2, and 
Q{b = ±90°) = 2.44 X IQ-^gcm-^. In order to 
compare with the EGRET results of Hunter et al. 
(1997), we have also evaluated the average of Q 
over the inner Galactic disk: {Q{\1\ < 60°, |6| < 
10°)) =3.28 X 10-2gcm-2. 

3. Discussion 

3.1. Comparison with data 

Several physical processes contribute to the ob- 
served diffuse gamma-ray intensity. In ordc;r of de- 
creasing fractional contribution these are thought 
to be pion production and bremsstrahlung from 
cosmic-ray interactions with diffuse Galactic gas, 
inverse Compton emission from cosmic-ray elec- 
trons interacting with ambient photons, and an 
isotropic background, which is presumably extra- 
galactic and due to many faint, discrete sources. 
The sum of these contributions offers a good model 
for the diffuse emission which is actually observed 
in the 100 MeV-1 GeV band (Bloemen 1989; 
Hunter et al. 1997); however, the E > 1 GeV 
emission is poorly modeled (Hunter ct al. 1997) 
with the prediction (Bertsch et al. 1993) amount- 
ing to only ~ 2/3 of the observed intensity. 

Some authors (Gilmore 1994; Salati et al. 
1996) have used the low-energy {E < 1 GeV) data 
to argue that the agreement between model and 
data allows no room for any significant unmod- 
eled emission and that these data therefore place 
tight constraints on any baryonic contribution to 
the dark matter halo of the Galaxy. This line of 
argument is clearly suspect because of the failure 
of the same model to account for the high-energy 
{E > 1 GeV) data. However, even if we ignore 
the high-energy data the calculations presented in 
§2.3 demonstrate that the gamma-ray constraints 
on high column density gas clouds in the dark 
halo are much weaker than those on low column 
density gas because the emissivity of the latter is 
much greater; consequently, we revisit the pub- 
lished constraints in the following section (§3.2), 
returning to the question of the high-energy data 
in §3.3. 

3.2. Published constraints applied to high 
column-density gas 

Gilmore (1994) argued that the known (diffuse) 



gas accounts for essentially all of the gamma-ray 
intensity in the E > 100 MeV band observed 
by the COS B satellite (Bloemen 1989), and 
he suggested that any contribution from bary- 
onic clouds should amount to no more than 
lO-^phcm-^s-^sr-i {E > 100 MeV) at high 
Galactic latitudes. Toward the Galactic poles 
we found (§3) Q ~ 2.5 x 10~^gcm~^, implying 
that the average emissivity of the material in the 
dark halo should he £ < 5.0 x 10~^phs~^g~^ for 
E > 100 MeV. From Table 1 (or figures 6 and 
7) we see that this requirement is not met for the 
models which are in the thin material limit but 
that any model with S > lOOgcm"^ is accept- 
able. In other words the constraint imposed by 
Gilmore (1994) permits the Galaxy's dark halo to 
be entirely baryonic, provided the individual com- 
ponents have a column-density S > lOOgcm"^. 

Salati ct al. (1996) pointed out that oberva- 
tions at low Galactic latitudes can provide tighter 
constraints on any baryonic dark halo, because 
of the higher intensity expected when looking 
through the cosmic-ray disk edge-on. This point 
is manifest in the large values oi Q{b — 0), which 
we computed in §2.3 (see also figure 8). Salati 
et al. (1996) argued that a suitable constraint 
on any baryonic component of the dark halo is 
that it should not contribute more than the uni- 
form background intensity component observed in 
any given field. The strongest constraint then 
comes from observations of the Ophiuchus region 
(Hunter et al. 1994), for which the requirement 
imposed by Salati et al. (1996) is an intensity 
contribution from the dark matter halo of Id < 
2.4 X 10-^phcm-=^s-^sr-i in the band 300 MeV 
<E< 500$ MeV. For this line of sight our calcu- 
lation yields Q{0,b = 15°) = 1.66 x lO'^gcm-^, 
and in the thin material limit, for which we find 
f ~ 2.4 x 10-2phs-ig-i (300 MeV < E < 500$ 
MeV; see Table 1), this corresponds to a predicted 
intensity of Id — 3.2 x 10~^ phcm~^s~^sr~^ in the 
300-500 MeV band, implying that at most 8% of 
the dark halo can be resident in low column den- 
sity clouds.^ However, our calculations extend to 

^This result is a factor of 2 larger than the correspond- 
ing model of Salati et al. (1996) (with L = 3 kpc and a 
spherical halo), a difference which is accounted for by their 
choice of a smaller core radius for the dark halo. Salati 
et al. (1996) chose a core radius of 3.5 kpc, which yields 
Q{0,b = 15°) = 3.0 X 10~'^gcm~'^, whereas we have em- 
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clouds of higher column-densities, and wc find that 
for S > 500gcm~^ the limit relaxes to the point 
where all of the dark halo is permitted to be in the 
form of high column density clouds. 

3.3. Revision of constraints based on the 
observed spectrum 

3.3.1. Galactic diffuse gamma-ray model 

The constraints discussed in the previous sub- 
section make reference only to the low-energy 
< 1 GeV) gamma-ray data. As noted earlier 
in this section, the observed intensity of the inner 
Galactic disk at high photon energies (> 1 GeV) is 
substantially greater than expected (Hunter et al. 
1997; Bertsch et al. 1993). Much effort has been 
expended on explaining this discrepancy, with 
most of the attention given to models in which 
the Galactic cosmic-ray electron and/or proton 
spectra differ from their locally measured values 
(Biisching, Pohl & Schlickeiser 2001; Aharonian 
& Atoyan 2000; Berezhko & Volk 2000; Strong, 
Moskalenko and Reimer 2000; Mori 1997). How- 
ever, none of these models offers a satisfactory ex- 
planation of the observed mean gamma-ray spec- 
trum of the Galactic disk, and consequently a suc- 
cessful match to the low energy data should not be 
taken to mean that the emission model is correct. 
In turn this suggests that the constraints formu- 
lated by Gilmore (1994) and Salati et al. (1996), 
on the basis of the low-energy data, may be too 
restrictive. 

The question then arises as to what constraints 
the gamma-ray data do in fact place on unniodelod 
emission, given the current state of understanding 
of the observed emission. 

The Galactic diffuse emission model used in 
Hunter et al. (1997) is based on Bertsch et al. 
(1993), and contains the following contributions: 

/ = /hi + /hii + + ^ic + /eg, (7) 

where /hi is the gamma-ray intensity contributed 
by cosmic-ray interactions with diffuse atomic hy- 
drogen, and similarly for the ionised and molecular 
components of the interstellar medium (HII and 
H2, respectively). The emissivity for these com- 
ponents (Bertsch et al. 1993) is, of course, com- 
puted in the thin material limit, and each compo- 

ployed a core radius of 6.2 kpc (see §2.3). 



nent therefore has the same spectral shape, differ- 
ing only in normalisation. Here he is the gamma- 
ray flux by inverse Compton emission (cosmic-ray 
electrons up scattering low-energy photons), and 
Ieg is the isotropic extragalactic background flux. 

Although the atomic and ionised hydrogen com- 
ponents can be observed directly via their line 
emission, and are thus well constrained, this is not 
true for the molecular component. The molec- 
ular hydrogen column is assumed to be propor- 
tional to the CO emission line strength, as mea- 
sured by the Columbia CO survey (Dame et al 
1987), for example, but the constant of propor- 
tionality (usually denoted by X) is unknown and 
one is forced to determine its value by fitting to 
the gamma-ray data. Although the uncertainty in 
the best-fit determination of {X) (averaged over 
the whole sky) is small, the systematic uncertain- 
ties are acknowledged to be much larger, "at least 
10%-15%" (Hunter et al. 1997). In turn, this 
estimate of the uncertainty is small in compari- 
son with the differences among the various values 
of X which have been reported in the literature 
(see the review by Bloemen (1989)) and the likely 
range of variation in X within the Galaxy (Hunter 
et al. 1997). Although these are substantial un- 
certainties, molecular hydrogen contributes only a 
fraction of the total observed gamma-ray intensity 
roughly 20% of the local surface density of gas in 
the Galactic disk is thought to be in molecular 
form (Wouterloot et al 1990) so the implied frac- 
tional uncertainty in the total Galactic emission is 
perhaps as small as 3%. Moreover, uncertainty in 
X affects only the normalisation and not the spec- 
tral shape of the predicted emission from diffuse 
molecular hydrogen, so the observed high-energy 
excess cannot be explained in this way even if X 
could assume an arbitrarily large value. 

The other main contributions to uncertainty in 
the diffuse model are associated with (1) unmod- 
eled spatial variations in the cosmic-ray spectral 
energy density and (2) unmodeled spatial vari- 
ations in the low-energy photon spectral energy 
density; the former affects all of the Galactic con- 
tributions to /, whereas the latter affects only the 
inverse Compton component. These uncertain- 
ties affect both the normalisation and the spectral 
shape of the predicted gamma-ray emission; how- 
ever, the uncertainties are difficult to quantify. 
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3.3.2. Constraints based on the gamma-ray data 

For OUT purposes it is not actually necessary 
to quantify the uncertainties on the model input 
parameters; it suffices to use the discrepancy be- 
tween model and data as a measure of the uncer- 
tainty in our understanding of the observed emis- 
sion. In turn this measure determines the con- 
straints which we can apply to any putative un- 
modeled emission, such as the contribution from 
dense gas which we are concerned with here. At 
photon energies E > \ GeV the fractional dis- 
crepancy is roughly 60% (Hunter et al. 1997), in 
the sense that the observed emission is 1.6 times 
larger than the model, and we henceforth adopt 
0.6/1.6 ~ 40% of the total observed intensity as 
our estimate of the unmodeled emission. Although 
this estimate is derived from data at high energies, 
the effects of the various contributing processes are 
all very widely spread, and the estimate therefore 
applies independent of photon energy. The con- 
straints appropriate to high/low Galactic latitudes 
can now be re-evaluated. 

At high Galactic latitudes the observed in- 
tensity is / ~ 1.5 X 10~^ phcm~^s~^sr~^ for 
E > 100 MeV (Kniffen et al. 1996), im- 
plying that any unmodeled emission should be 
< 6 X 10~^ phcm^^s^^sr"^ in this band. This 
result is actually slightly stricter than the crite- 
rion used by Gilmore (1994) and thus leads us to 
tighten our high-latitude constraints, relative to 
those quoted in §3.2: the observed high-latitude 
gamma-ray intensity constrains the amount of low 
column- density gas to < 20% of the total den- 
sity of the Galactic dark halo, with this fraction 
rising to 100% for gas clouds of column density 
S > 200gcm-2. 

At low Galactic latitudes we can make use 
of the mean intensity of the inner Galactic 
disk, which has been accurately determined by 
Hunter et al. (1997). For example at 1 GeV 
the mean intensity < 60°, |6| < 10°) is 
(/) ~ 3 X 10"*phcm"2s"isr"iMeV"\ and our 
calculation of (Q) for this region yields (§2.3) 
3.28 X lO^^gcm^^, implying that the emissiv- 
ity of the Galactic dark halo material must be, 
on average, f < 4.6 x 10^® phs^^g^^MeV^^. By 
comparison, the actual emissivity of low column- 
density gas is computed to be (§2.3, table 2) 
£{1 GeV) ~ 1.4 x 10-5phs-ig-^MeV-\ imply- 



ing that < 30% of the Galaxy's dark halo may be 
comprised of low column density gas. For higher 
column densities the emissivity falls, and table 2 
shows that for S = 100 gcm^^ the emissivity is 
only 5.5 x 10~^ phs~^g~^MeV~^. The gamma-ray 
data on the inner Galactic disk thus indicate all of 
the Galaxy's dark halo to be made of dense clouds 
of column-density S > lOOgcm"^. 

3.4. Comment on the gas content of the 

galactic disk 

The constraint we have just given is based on 
the mean spectrum of the inner Galactic disk, in 
contrast to those given by Salati et al. (1996) 
who employed limits based on the angular struc- 
ture of the observed gamma-ray intensity. Specif- 
ically, Salati et al. (1996) required that the puta- 
tive contribution of emission from a baryonic com- 
ponent of the Galaxy's dark halo be less than that 
of the isotropic component of the observed inten- 
sity; this procedure seems to us to be less reli- 
able than the procedure we have employed, for two 
reasons. First, even if the dark halo were spher- 
ically symmetric the emission from any baryonic 
component would not be, both because our point 
of observation is quite distant from the centre of 
the Galaxy and because the resulting gamma-ray 
emission is strongly dependent on the Galactic 
cosmic-ray distribution, which in turn is strongly 
concentrated in the disk of the Galaxy. Indeed the 
cosmic-ray distribution appears to correlate with 
the distribution of interstellar matter (Hunter et 
al. 1997), thus complicating the interpretation of 
the observed gamma-ray intensity distribution. In 
particular this coupling leads to gamma-ray emis- 
sion from a baryonic dark halo being correlated 
with the diffuse gas column density, even if the 
dense gas is uncorrclated with the diffuse gas. Sec- 
ond, on any given line of sight, such as the Ophi- 
uchus field considered by Salati et al. (1996), a 
highly structured dark matter halo might exhibit, 
by chance, a low dark matter column density. The 
chances of this are good if, as in the case of Salati 
et al. (1996), the field is specifically chosen to have 
a low "background" intensity. 

4. Summary 

The gamma-ray spectra arising from cosmic-ray 
interactions with gas clouds of various column- 
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densities have been calculated using a Monte Carlo 
event simulator, GEANT4. Our calculations re- 
produce the analytic result in the low column- 
density limit, where only single particle interac- 
tions need to be considered, but exhibit signifi- 
cant differences for clouds of column-density E > 
10^ gcm~^ where the emissivity declines substan- 
tially for photon energies E > 100 MeV. The low 
emissivity of dense gas means that the baryonic 
content of the Galaxy's dark halo is not so tightly 
constrained by the gamma-ray data as had previ- 
ously been thought. For S > 200 gcm~^ we find 
that the existing gamma-ray data, taken in isola- 
tion, do not exclude purely baryonic models for 
the Galactic dark halo. 

We are grateful to the referee for constructive 
comments that have improved this paper. 
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Kinetic energy, Tp (GeV) 



Fig. 1. — Adopted cosmic-ray proton spectrum 

in the Galaxy. The model curves are taken from 
Mori (1997) and Shikaze (2001). Also plotted are 
direct observations of the local spectrum (Ryan et 
al. 1972; Smith et al. 1973; Webber et al. 1987; 
Kawamura et al. 1989; Seo et al. 1991; Ivanenko 
et al. 1993; Asakimori et al. 1998; Boezio et al. 
1999; Sanuki et al. 2000), but note that these are 
strongly affected by solar modulation. 




Moskalenko and Strong 98 {w/ reacc) 
— Casadei and Bindi 03 



10'^ 10''' lo" 10'' 10^ 10^ 10'' 

Kinetic energy, T, (GeV) 

Fig. 2. — Adopted cosmic-ray electron spectrum 
in the Galaxy. The model curves are taken 
from Skibo and Ramaty (1993), Moskalenko and 
Strong (1998) and Casadei and Bindi (2003). 
(The two lines for Moskalenko and Strong (1998) 
show models with and without reacceleration of 
the cosmic rays.) Also plotted are direct observa- 
tions of the local electron spectrum (Tang 1984; 
Golden et al. 1984, 1994; Barwick et al. 1998; 
Kobayashi et al. 1999; Boezio et al. 2000; Torii 
et al. 2001), but note that these are affected by 
solar modulation. 
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lO'' 10^ 10^ lO"* 10^ 10^ 10^ 
Gamma-ray energy (IWeV) 

Fig. 3. — Gamma-ray emissivity, , for a molec- 
ular hydrogen cloud, of radius 1 AU and column 
density 100 gcm~^, irradiated by cosmic-ray pro- 
tons; vertical bars represent Monte Carlo statisti- 
cal errors. The various contributing emission pro- 
cesses are also shown separately. 
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Fig. 4. — Gamma-ray emissivities for dense clouds 
irradiated by cosmic-ray protons, with a spec- 
trum appropriate to the solar neighbourhood. The 
results have been multiplied by a nuclear en- 
hancement factor of 1.52 (Mori 1997) to account 
for the presence of heavier nuclei in the cosmic- 
rays. The different curves are for various different 
column-densities: 10~^gcm~^, thick dot-dashed 
curve; 10~^gcm~^, thick dashed curve; lgcm~^, 
thick dotted curve; lOgcm"^, thick solid curve; 
lO^gcm"^, thin dot-dashed curve; lO^gcm"^, 
thin dashed curve; 10'*gcm~^, thin dotted curve; 
lO^gcm"^, thin solid curve. Also plotted is the 
emissivity from Mori (1997) , which corresponds to 
the "thin material" limit (filled circles); this limit 
offers a good approximation for E < lOgcm"^. 
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Fig. 5. — Gamma-ray emissivities for dense clouds 
irradiated by cosmic-ray electrons. The different 
curves are for different column-densities, as per 
figure 4. 
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Fig. 6. — Integrated gamma-ray emissivities for 
dense clouds irradiated by cosmic-ray hadrons. 
Solid curves are those assuming a cosmic ray pro- 
ton flux in Mori (1997) and dotted ones assuming 
that in Shikaze (2001). Here we multiplied by a 
nuclear enhancement factor of 1.52 (Mori 1997) 
to take account of heavy nuclei in the cosmic rays. 
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Column density (g om'^) 

Fig. 7. — Integrated gamma-ray emissivities for 
dense clouds irradiated by cosmic-ray electrons. 
Solid curves arc those assuming a cosmic ray elec- 
tron flux in Skibo and Ramaty (1993), and dashed 
(dotted) ones assuming that in Moskalenko and 
Strong (1998) with (without) re-acceleration. 




Fig. 8. — Adopted model dark matter column- 
density distribution over the sky, Q{l,b), weighted 
by cosmic-ray energy density and displayed in 
Galactic coordinates. The quantity Q is de- 
fined in equation 4 and is computed using 
the dark matter density model given in equa- 
tion 5 together with the cosmic-ray distribution 
model given in equation 6. Contour levels are 
(4, 8, 12, 16, 20, 24, 32, 40, 48, 56, 64) x IQ-^ g cm-^; 
dashed lines are constant longitude, spaced at 
60° intervals; dotted lines are constant latitude, 
spaced at 30° intervals . 
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